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Abstract 



This paper reports an updated measurement of the Standard Model CP 
violation parameter sin 2/3 using the CDF Detector at Fermilab. The entire 
Run I data sample of 110 pb~^ of proton antiproton collisions at =1.8 TeV 
is used to identify a signal sample of ~400 B J/^Kg events, where J/V' — 
fi'^H" and Kg tt'^tt'. The flavor of the neutral B meson is identified at the 
time of production by combining information from three tagging algorithms: a 
same-side tag, a jet-charge tag, and a soft-lepton tag. A maximum likelihood 
fitting method is used to determine sin 2(3 = 0.79+°;44(stat-fsyst). This value 
of sin 2/3 is consistent with the Standard Model prediction, based upon existing 
measurements, of a large positive CP-violating asymmetry in this decay mode. 

PACS numbers: 12.15.Hh, 13.20.He, 14.40.Nd 
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I. INTRODUCTION 



The first observation of a violation of charge-conjugation parity (CP) invariance was 
in the neutral kaon system in 1964 [|I|. To date, violation of CP symmetry has not been 
directly observed in any other system. The study of CP violation in the B system is an ideal 
place to test the predictions of the Standard Model The decays of neutral B mesons 

into CP eigenstates are of great interest, in particular the CP-odd state, B —> J/ipKg [^0. 
The decay B —>■ J/tpK^ is a popular mode in which to observe a CP- violating asymmetry 
because it has a distinct experimental signature and is known theoretically to be free of large 
hadronic uncertainties 0. Furthermore, the contribution to the asymmetry due to penguin 
diagrams, which is difficult to calculate, is negligible because the penguin contribution is 
small and the tree level and penguin diagrams contribute with the same weak phase 0. 
Previous work searching for a CP- violating asymmetry in the decay B J/ipKg has been 
presented by the OPAL Collaboration An initial study on the measurement of sin 2/3 
by the CDF Collaboration is given in Ref. |T^. The result reported here incorporates and 



supersedes Ref. [|T0|. This paper reports a measurement of sin 2/3 that is the best direct 
indication of a CP- violating asymmetry in the neutral B meson system. 

Within the framework of the Standard Model, CP nonconservation arises through a 



non-trivial phase in the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix |]TT 
The CKM matrix V is the unitary matrix that transforms the mass eigenstates into the 
weak eigenstates: 



V 



' Vud Vus Vub 
Vcd Vcs Vcb 

V Vtd Vt, Vtb 
1 - 

I 

p-iri) 





The second matrix is a useful phenomenological parameterization of the quark mixing matrix 
suggested by Wolfenstein |12[, in which A is the sine of the Cabibbo angle. The condition 



of unitarity, V'^V = 1, yields several relations, the most important of which is a relation 
between the first and third columns of the matrix, given by: 

v:bVud + v:,v^ + v:,Vt, = o. 

This relation, after division by V^Vcd, is displayed graphically in Fig. |I] as a triangle in 
the complex {p-rj) plane, and is known as the unitarity triangle |13|. CP violation in the 



Standard Model manifests itself as a nonzero value of r], the height of the triangle. 

CP nonconservation is expected to manifest itself in the B^ system p[ as an asymmetry 
in particle decay rate versus antiparticle decay rate to a particular final state: 

^CP — =n 

N{B ^ J/^K'^s) + ^(5° ^ J/i^K'^s) 

where J/ipK^) is the number of mesons decaying to J/ipKg that were produced as 

B and N{B^ J/ipKg) is the number of mesons decaying to J/ipKg that were produced 
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as 5° 0. It should be noted that the definition of Acp is the negative of that in Refs. 
and H. 

In the Standard Model, the CP asymmetry in this decay mode is proportional to sin 2/?: 
Acp{t) = sin 2/? sin(Am(it), where f3 is the angle of the unitarity triangle shown in Fig. |l|, 
t is the proper decay time of the meson and Arrid is the mass difference between the 
heavy and light mass eigenstates. In a hadron collider, BB pairs are produced as two 
incoherent meson states. Consequently, the asymmetry can be measured as either a time- 
dependent or time-integrated quantity. The time- dependent analysis is however statistically 
more powerful. In this paper, we take advantage of this fact and employ a sample of events 
that have a broad range of time resolutions. 

It is possible to combine information from several measurements to indirectly constrain 
the allowed range of sin 2/5. Based on global fits to these measurements, it is found that 
the Standard Model prefers a large positive value of sin 2/5 and that the fits are in good 
agreement with each other ||l^^l7|. One recent global fit finds sin 2/5 = 0.75 ± 0.09 [l7 



However, the sign of the expected asymmetry depends on the sign of the product of Bb and 
Bk, which are the ratios between the short distance contributions to BB and KK mixing 



respectively and their values in the vacuum insertion approximation [|TI 

To measure this asymmetry, the fiavor of the B meson (whether it is a i?*^ or a 5") must 
be identified (tagged) at the time of production. The effectiveness of a tagging algorithm 
depends on both the efficiency for assigning a flavor tag and the probability that the flavor 
tag is correct. The true asymmetry is "diluted" by misidentifying a B^ meson as a i? meson 
or vice versa. We deflne the tagging dilution a.s D = {Nji — N^) / {Nji + N^) , where Nb,{Nw) 
is the number of right (wrong) tags. The observed asymmetry, given by A'^p = DAcp, is 
reduced in magnitude by this dilution parameter. As can be seen from the relation above, 
maximal sensitivity to the asymmetry is achieved when the dilution factor is large. The 
statistical uncertainty on sin 2/5 is inversely proportional to velP, where the efficiency e is 
the fraction of events that are tagged. This analysis combines three tagging algorithms in 
order to minimize the statistical uncertainty of the measurement. 



A. The CDF detector 

The CDF detector is described in detail elsewhere |]TU|,PD[- The CDF detector systems 



that are relevant for this analysis are: (i) a silicon vertex detector (SVX) (ii) a time 
projection chamber (VTX), (iii) a central tracking chamber (CTC), (iv) electromagnetic and 
hadronic calorimeters, (v) a preshower detector (CPR, central preradiator) , (vi) a shower 
maximum detector (CES, central electron strip chamber), and (vii) a muon system. The 
CDF coordinate system has the 2;-axis pointing along the proton momentum, with the x-axis 
located in the horizontal plane of the Tevatron storage ring, pointing radially outward, so 
that the y-axis points up. 

The SVX consists of four layers of silicon axial-strip detectors located between radii of 
2.9 and 7.9 cm and extending ± 25 cm in z from the center of the detector. The geometrical 
acceptance of the SVX is ~ 60% because the pp interactions are distributed with a Gaussian 
profile along the beam axis with a standard deviation of ~ 30 cm, which is large relative to 
the length of the detector. The SVX is surrounded by the VTX, which is used to determine 
the z coordinate of the pp interaction (the primary vertex). Momenta of charged particles 
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are measured in three dimensions using the CTC, an 84-layer drift chamber that covers the 
pseudorapidity interval \r]\ <1.1, where t] = — ln[tan(^/2)], and the angle 6 is measured 
from the z-axis. The SVX, VTX, and CTC are immersed in a 1.4 T solenoidal magnetic 
field. The momentum transverse to the beamline (Pt) of a charged particle is determined 
using the SVX and CTC detectors. The combined CTC/SVX Pt resolution is 6Pt/Pt = 
[(0.001 c/GeV ■ Pt)^ + (0.0066)^] 2. The typical uncertainty on the B meson decay distance 
is about 60 /xm. The CTC also provides measurements of the energy loss per unit distance, 
dE/dx, of a charged particle. 

The central and endwall calorimeters are arranged in projective towers and cover the 
central region |?7| <1.05. In the central electromagnetic calorimeter, proportional chambers 
(CES), are embedded near shower maximum for position measurements. The CPR is located 
on the inner face of the central calorimeter and consists of proportional chambers. The muon 
system consists of three different subsystems each containing four layers of drift chambers. 
The central muon chambers, located behind ~ 5 absorption lengths of calorimeter, cover 85% 
of the azimuthal angle (j) in the range \ri\ < 0.6. Gaps in are filled in part by the central 
muon upgrade chambers with total coverage in of 80% and |?7| < 0.6. These chambers 
are located behind a total of ~ 8 absorption lengths. Finally, the central extension muon 
chambers provide 67% coverage in (p for the region 0.6 < |?7| < 1.0 behind a total of ~ 6 
absorption lengths. 

Muons, used to reconstruct the J/ip meson and by the soft lepton tagging algorithm 
(SLT), are identified by combining a muon track segment with a CTC track. SVX infor- 
mation is used when available. Electrons, which are used by the SLT, are identified by 
combining a CTC track with information from the central calorimeters, the central strip 
chambers, dE/dx, and the CPR detectors. 

Dimuon events are collected using a three-level trigger. The first-level trigger system 
requires two charged track segments in the muon chambers. The second level trigger requires 
a CTC track, with Pt greater than ~ 2 GeV/ c, to match a muon chamber track segment. The 
third level, implemented with online track reconstruction software, requires two oppositely 
charged CTC tracks to match muon track segments and a dimuon invariant mass between 
2.8 and 3.4 GeV/c^. Approximately two thirds of all J/ip ^ events recorded enter 

on a dedicated J/ip trigger, where the two reconstructed muons are from the J/ip- This 
fraction is consistent with expectations. The majority of the remaining events, referred to 
as "volunteers", enter the sample through: a single inclusive muon trigger caused by one 
of the two muons from the J/ip decay, or, through a dimuon trigger where one of the two 
trigger muons was from the J/iJj and the second "trigger muon" is a fake muon, primarily 
due to punch-through. 



B. Overview of the analysis 

This analysis builds on the work of several previous analyses using the various B enriched 
data sets recorded by the CDF detector. The B J/ipK^ decay mode is reconstructed 
in a manner similar to the CDF measurements of the branching ratio ||2^j23|] and the B 
lifetime |^^. The three tagging algorithms are then applied to the B J/ipKg sample 



and the observed asymmetry, given by A^p = DAqp, is then determined. In order to 
extract a value of sin 2/? from the observed asymmetry, tagging dilution parameters are 
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required for the three tagging algorithms. These dilution parameters are determined from 
an analysis of the calibration samples. In particular, the same-side tagging (SST) dilutions 
are determined from a combination of results from Ref. and measurements on a sample 
of ~1000 -B^ J/ipK^ decays. The jet-charge tag algorithm (JETQ) and soft-lepton tag 
algorithm (SLT) dilutions are determined from the J/ipK^ sample and ~ 40, 000 

inclusive B J /ipX events. The dilutions and efficiencies are then combined for each event 
and a maximum likelihood fitting procedure is used to extract the result for sin 2/?. The fit 
includes the possibility that the tagging dilutions and efficiencies have inherent asymmetries. 
In addition, the backgrounds, divided into prompt and long-lived categories, are also allowed 
to have an asymmetry. In the end, these possible asymmetries are found not to be significant. 

Each fiavor tagging method, SST, SLT and JETQ, has been previously verified in a B^- 
B mixing analysis. Our previously published measurement of sin 2(3 used the B^-B mixing 
analysis of Ref. ||25| to establish the viability of the SST method [^]. Here we report work 
that uses the same algorithm for events where the two muons are contained within the SVX 
detector acceptance and uses a modified version of the algorithm for events with less precise 
fiight path information, i.e. events not fully contained within the SVX detector acceptance. 

The two additional tagging algorithms used are based on the B^-^ mixing analysis of 
Ref. These mixing analyses use decays of B mesons with higher (~ a factor of 

two higher) than the B mesons in this analysis. This is due to the lower trigger threshold 
for J/ip ^ fi^fi~ than for the inclusive lepton triggers used to select the mixing analyses 



samples. The SLT algorithm is similar to that in Ref. [27], except the lepton Pt threshold has 



been lowered to increase the efficiency of tagging lower Pt B mesons. The JETQ algorithm 
is also similar to the algorithm used in the mixing analysis |27] except the acceptance cone 
defining the jet has been enlarged and impact parameter weighting of tracks has been added 
to reduce the fraction of incorrectly tagged events. 



II. SAMPLE SELECTION 

Four event samples, B J/ipKg, B^ J/ipK^, inclusive B J/ipX decays, and 
an inclusive lepton sample |^ are used in the determination of sin 2/3. The B mesons are 
reconstructed using the decay modes J/V' — > /^"""/i" and Kg — > tt+tt". The B — >• J/ipK^ 
candidates form the signal sample, the B^ J/ipK^ sample is used to determine the 
tagging dilutions, and the inclusive J/ip decays are used to constrain ratios of efficiencies. 
The inclusive lepton sample was used in Refs. [p!0| , p5| in the determination of the SST 
dilution. 

The selection criteria are largely the same as in Ref. The criteria for the B J/ipKg 
sample provide an optimal value of the ratio 5'^/(5'-|-iVbck), where S is the number of signal 
events and A^bck is the number of background events within three standard deviations of 
the B mass. The square root of this ratio enters into the uncertainty on the measurement 
of sin 2/5. The J/ip is identified by selecting two oppositely charged muon candidates, each 
with Pt > 1.4 GeV/c. Additional selection criteria are applied to ensure good matching 
between the CTC track and the muon chamber track segment. A J/ip candidate is defined 
as a pair within ±5cr of the world average mass of 3.097 GeV/c^ [Q, where a is the 

mass uncertainty calculated for each event. 



8 



The Kg candidates are found by matching pairs of oppositely charged tracks, assumed 
to be pions. The candidates are required to travel a significant distance L^y > 5a l, 
and to have Pt > 700 MeV/c in order to improve the signal-to-background ratio. The 
quantity L^y = X ■ P^- is the 2-D flight distance, where X is the vector pointing from 
the production vertex to the decay vertex, and aL is the measurement uncertainty on L^y. 
This flight distance is used to calculate the proper decay time t = L^yMo/ Pt, where Mq 
is the world average 5° mass of 5.2792 GeV/c^ §. In about 15% of the /sTg decays, SVX 
information is available for one or both tracks. When the decay vertex location in the 
radial direction is found to lie beyond the second layer of the SVX detector, the SVX 
information is not used. The J/ip and Kg candidates are combined into a four particle fit 
to the hypothesis B — > JjtpKg and the /i"*"/!" and tt+tt^ are constrained to the appropriate 
masses and separate decay vertices. The K^ and B are constrained to point back to their 
points of origin. In order to further improve the signal-to-background ratio, B candidates 
are accepted for Pt{B) > 4.5 GeV/c and fit quality criteria are applied to the J/ip and B 
candidates. 

The data are divided into two samples, one called the SVX sample, the other the non- 
SVX sample. The SVX sample requires both muon candidates to have at least three out of 
four possible hits that are well measured by the silicon vertex detector. This is the sample 
of B candidates with precise decay length information and is similar to the sample that was 
used in the previously published CDF sin 2/5 analysis. The non-SVX sample is the subset of 
events in which one or both muon candidates are not measured in the silicon vertex detector. 
About 30% of the events in this sample have one muon candidate track with high quality 
SVX information. Events of this type lie mostly at the boundaries of the SVX detector. 

We define a normahzed mass M^v = (^/^^tttt — Afo)/c"fit, where m^^^^ is the four-track 
mass coming from the vertex and mass-constrained fit of the B candidate. The uncertainty, 
(jfit, is from the fit, typically ~ 10 MeV/c^. The normalized mass distribution is shown 
in Fig. ^ and contains 4156 entries, from which we observe 395 ± 31 signal events with a 
signal-to-noise ratio of 0.7. The SVX sample contains 202 ± 18 events (signal-to- noise ratio 
of 0.9) and the non-SVX sample contains 193 ± 26 events (signal-to-noise ratio of 0.5) as 
shown in Fig. ^. The event yields reported here come from the full unbinned likelihood fit 
which will be described in detail later. 

The criteria used to select the B^ —* J/ipK^ decays are the same as described for 
B J/ipKg decays except for the K"^ selection. Since the CDF detector has limited 
particle identification separation power at high using the dE/dx system, candidate kaons 
are defined as any track with P^ > 2 GeV/c. The fi^fi^K^ mass distribution is shown in 
Fig. ^ and the number of J/ipK^ candidates is 998 ± 51. 

The inclusive J/ip ^ sample is a superset from which the B J/ipKg and 

P^ — i> J/tpK^ samples are derived. The inclusive sample is ~80% prompt J/ip from direct 
cc production. In order to enrich the sample in P — J/ ipX decays, both muons are required 
to have good SVX information and the J/ ip 2-D travel distance must be > 200 ixm from the 
beamline. This results in a sample of about 40,000 P J/ipX decays. 



9 



III. TAGGING ALGORITHMS 



Three tagging algorithms are used, two opposite-side tag algorithms and one same-side 
tag (SST) algorithm. The idea behind the SST algorithm [^] exploits the local correla- 



tion between the B meson and the charge of a nearby track to tag the flavor of the B 



meson. We employ the SST algorithm described in detail in Ref. [10,25|. We consider all 



charged tracks that pass through all stereo layers of the CTC and within a cone of radius 

AR = y/Arf~+~A(f^ < 0.7 centered along the B meson direction. Candidate tracks must 

be consistent with originating from the primary vertex and have a. Pt > 400 MeV/c. If 

more than one candidate is found, the track with the smallest P|?^ is chosen, where P|?^ is 

the track momentum transverse to the momentum sum of the track and the B meson. A 

— 

tagging track with negative charge indicates a B meson, while a positive track indicates a 
B^ meson. 

The performance of the SST algorithm could depend on the availability of precise vertex 
information. When using the SVX sample, the SST algorithm of Ref. |Ty] and tagging dilu- 
tion parameter D = (16.6 ± 2.2)% is used. This dilution result is obtained by extrapolating 
the value obtained in the mixing analysis in Ref. to the lower Pt of the B J/ipKg 
sample. When using the non-SVX sample, the SST algorithm is modified slightly by drop- 
ping the SVX information for all candidate tagging tracks and adjusting the track selection 
criteria in order to increase the geometrical acceptance. A dilution scale factor Jd, defined 
by -D non-SVX = /d-Dsvx, is derived from the B^ J/ipK^ sample. This relates the SVX 
sample SST algorithm performance to that of the non-SVX sample SST algorithm. To mea- 
sure this quantity, we compare the tagging track using SVX information to the track we 
obtain when all SVX information is ignored. This provides a measure of the effectiveness of 
the SVX information. We find a value of fc = (1.05 ± 0.17), apply it to the measured SST 
dilution for SVX tracks, and obtain D = (17.4 ± 3.6)%. 

Opposite-side tagging refers to the identification of the flavor of the "opposite" B in 
the event at the time of production. As mentioned earlier, two algorithms are employed: 
soft-lepton tag (SLT) and jet-charge tag (JETQ) algorithms. 

The SLT algorithm is described in detail in Ref. [^]. The SLT algorithm associates 
the charge of the lepton (electron or muon) with the flavor of the parent S-meson, which 
in turn is anticorrelated with the produced flavor of the i?-meson that decays to J/ijjKg. 
These leptons are considered "soft" because their momenta are on average considerably 
lower than the high momentum leptons from W boson, Z boson, and top quark decays. A 
soft muon tag is defined as a charged track reconstructed in the CTC (CTC track) with 
Pt > 2 GeV/c that has been matched to a track segment in a muon system. A soft electron 
tag is defined as a CTC track with Pt > 1 GeV/c that has been successfully extrapolated 
into the calorimeters, CPR and CES detectors and passed selection criteria. In particular, 
the CPR and CES position information is required to match with the CTC track and the 
shower profiles must be consistent with an electron. In addition, the electron candidate 
CTC track must have a dE/dx deposition consistent with an electron. Photon conversions 
are explicitly rejected. A dilution of -D = (62.5 ± 14.6)% is obtained by applying the SLT 
algorithm to the -B^ — > J/ipK^ sample. 

If a soft lepton is not found, we try to identify a jet produced by the opposite B. We 
calculate a quantity called the jet charge Qjet of this jet: 
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^j^* E.Pt.(2-(T,),) ' 

where qi and P^j are the charge and transverse momentum of the i^^ track in the jet with 
Pt > 750MeV/c. The quantity Tp is the probabihty that track i originated from the pp 
interaction point. The quantity (2 — Tp) is constructed such that a displaced (prompt) track 
has the value Tp ~0 (1), and the quantity (2 — Tp) is ~2 (1). Tracks that arise from B decays 
are displaced from the primary vertex and give a probability distribution Tp peaked near 
zero, lending larger weight to the sum. For tracks that emanate from the primary vertex, 
Tp is a flat distribution between and 1, giving less weight to the jet charge quantity. For 
6-quark jets, the sign of the jet charge is on average the same as the sign of the 6-quark 
that produced the jet, so the sign of the jet charge may be used to identify the flavor at 
production of the B hadron which decayed to J/ipKg. This algorithm is conceptually similar 



to that used in Ref. |23] except that jet clustering and weighting factors are optimized for this 
sample. This optimization was performed by maximizing eD^ on a sample of B"^ J/ipK^ 
events generated by a Monte Carlo program. 

Jets are found with charged particles instead of the more commonly used calorimeter 
clusters. The algorithm is optimized using Monte Carlo generated data. All tracks in an 
event with Pt > 1.75 GeV/c are identified as seed tracks. For pairs of seed tracks, the 
quantity Yij = 2EiEj{l — cos 9ij) is calculated, where Ei, Ej are the energies and 9ij is the 
angle between the ith and jth seed tracks. Seed tracks are combined in pairs as long as Yij, 
the JADE distance measure, is less than 24 GeV^. After mergings, each set of seed tracks 
defines a jet. The remaining tracks [Pt < 1.75 GeV/c) are combined with the jet that 
minimizes the distance measure provided that Yij < 24 GeV^. Any tracks unassociated with 



a track-group are discarded. This is a modified version of the JADE clustering algorithm pS 

Tracks within a cone of AR < 0.7 with respect to the B — ^ J/ipKg direction are excluded 
from clustering to avoid overlap with the SST candidate tracks. The B meson decay products 
(/i"'",/i~,7r"'" and vr") are also explicitly excluded from the track-group. A jet can consist of 
a single track with Pt > 1.75GeV/c. If multiple jets are found, we choose the one that is 
most likely a 5 jet, based on an algorithm that uses the track impact parameter information 
first, if available, and then the jet Pt- The momentum and impact parameter weighted 
charge, Qjct, is calculated for the jet and normalized such that |Qjet| < 1- Only tracks with 
Pt > 0.750 GeV/c are used to weight the charge. The parameter Qjct > 0.2 selects the b 
quark decays and Qj^t < —0.2 selects the b quark decays. The value |Qjct| <0.2 is considered 
untagged. A dilution of D = (23.5 ± 6.9)% is found by applying the JETQ algorithm to the 
B^ —>■ J/tpK^ sample. 

We use a sample of 998 ± 51 B"^ J/ipK"^ decays to determine the tagging dilutions 
for the opposite-side algorithms. Using both real data and simulated data, we have verified 
that D(i?^) is consistent with D{B^) for the opposite-side flavor tagging algorithms. At 
the Tevatron, the strong interaction creates bb pairs at a production energy sufficiently high 
that the fragmentation processes that create the B mesons are largely uncorrelated. For 
example, the b quark could hadronize as a B~ meson, while independently, the b quark 
could hadronize as a B~^, B^ or B^ meson. These opposite side dilution numbers are valid 
for both the SVX and non-SVX samples. The tagging dilutions and efficiencies are presented 
in Table |. 

Each event has the opportunity to be tagged by two tag algorithms: one same-side and 
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one opposite-side. We followed the prescription outlined in Ref. in which the SLT tag 
is used if both the SLT and JETQ tags are available. This is done to avoid correlations 
between the two opposite side tagging algorithms. The result of the SLT algorithm is used 
because the dilution of the SLT algorithm is much larger than that of the JETQ algorithm. 
Given the low efficiency for lepton tags (6%) the potential overlap is small. As mentioned 
earlier, tracks eligible for the SST algorithm are excluded from the JETQ track list, thus 
ensuring these two algorithms are orthogonal. There is however an overlap between the SST 
and the SLT algorithms in which the lepton is used as the SST track. In order to use the 
dilution measured in Ref. we use the identical SST algorithm on the SVX sample, and 
therefore permit this overlap. We allow leptons in the cone to account for bb production 
from the higher-order gluon splitting process where the b — ^ iX decay is located nearby the 
fully reconstructed B —* J/ijjKg. This overlap occurs in three events in the signal region 
and the final result changes negligibly if these events are removed from the sample. 

Based upon the tagging efficiency of each individual tagging algorithm, we can calculate 
the expected fraction of events which will tagged by two, one or zero algorithms. We find 
the expected efficiency of each combination of tags {e.g. events tagged by both SST and 
SLT, events tagged by JETQ only, etc.) is consistent with estimates derived from a study 
of tagging efficiencies as applied to the B"^ J/ipK^ sample. Tag efficiencies are higher, 
typically by ~10%, in the trigger volunteer sample, except for the JETQ tagging algorithm, 
in which the efficiency increases by about 17%. These higher efficiencies are due to the 
increased average charged-track multiplicity of the trigger volunteer sample. Thus trigger 
samples that do not include volunteers, as planned for Run II, will have lower tagging 
efficiencies. It is found that ~ 80% of the events in the entire B J/ipKg sample are 
tagged by at least one tagging algorithm. 



A. Tag sign definition 

An event is tagged if it satisfies the criteria of any of the three tag algorithms. For all 
tag algorithms, the flavor tag refers to whether the candidate B J/ipKg was produced 
5° or B . The sign of all tag algorithms follow the convention established by the 



same-side tag algorithm discussed in Ref. |10|: The positive tag {+ tag) is deflned as the 



identiflcation of a 6-quark and therefore a B meson. The negative tag (— tag) is deflned as 
the identiflcation of a 6-quark and therefore a B^ meson. A null tag (or tag 0) means the 
criteria of the tag algorithms were not satisfled, and the flavor of the B is not identifled. A 
summary is provided in Table |l[ 



IV. DILUTIONS, EFFICIENCIES AND TAGGING ASYMMETRIES 

The dilutions and efficiencies described earlier need to be generalized in order to ac- 
commodate possible detector asymmetries in the analysis. For example, the CTC has a 
small (~ 1%) bias toward reconstructing more tracks of positive charge at low transverse 
momentum. This small bias is due to the tilted drift cell that is necessary to compensate 
for the Lorentz angle of the drift electrons, and a known asymmetry in background tracks 
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from beam pipe interactions. The formalism for measuring and correcting for these possible 
tagging asymmetries in this multi-tag analysis is provided below. 

For B mesons decaying to a CP eigenstate, the decay rate as a function of proper time 
t can be written as 

h±{t) = ——(1 ± Acpsin(Amdi)) 
It 

where hj^{t) is the decay rate for B's produced as type "+", h_{t) is the decay rate for B's 
produced as type "— ", and Acp = — sin 2/? is the asymmetry due to CP violation. Particle 
type "+" refers to a S ^ J/i/jKg decay and particle type "— " refers to a 5 — > J/ipK'g 
decay. 

To allow for an imperfect and (possibly) asymmetric tagging algorithm, the following 
definitions are used. For those B mesons of (produced) type +, a fraction will be actually 
tagged +, fraction will be tagged as — , and fraction eg will not be tagged, i.e. tag 0. 
Similarly, for those B mesons of (produced) type — , will be tagged — , fraction will be 
tagged as +, and fraction will be tagged as 0. Because, by definition, CR + e^^ + eo — 1 s-^id 
e]j + ej^ + eo = 1, there are 4 independent numbers that characterize a general asymmetric 
tagging algorithm. 

We define the efficiencies and dilutions for the general asymmetric tagging algorithm as 
e+ = (4 + evK)/2, e- = (cr + evF)/2, eo = {e^ + eg )/2 and 

n _ ^fl ~ n _ ~ n _ ~ 

^+ ~ + , - > ^- — - . + ■, ^0 — + _ ■ 

' ^R' ^0 ' ^0 

The observed decay rate as a function of time for events tagged as +, — or is given by 

g-t/r 

h+{t) = €+{1 + AcpD+sm{Amdt)), 



h^{t) = e_(l - Acp-D-sin(Amdi)), 



and 

^o(^) = eo{l + AcpDosm{Amdt)). 



Note that e+ + e_ + cq = 1 and e_|_D_|_ — e_D_ + eo^^o — 0) so there are four independent 
parameters remaining. For example, 

e_L>_ - e+D+ 



1 — e+ — e_ 
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A. Combining tags in an event 



Tagging information for each event is combined to reduce the uncertainty on the CP 
asymmetry. The tags are weighted for each event by the dilution of the individual tag 
algorithms. This procedure must also combine the efficiencies in a similar manner. The 
algorithm used to combine multiply-tagged events is as follows. Wc define the tags for 
two tagging algorithms as qi and q2 (each taking the values —1, 0, and 1), the individual 
dilutions as Di and D2, and the individual efficiencies as e^j and Cg^. We then define the 
dilution- weighted tags Vi = QiDi, the product of the tag and the dilution. We calculate the 
combined dilutions and efficiencies as 

= 1 + Dil)2 ^^^^^ ^ ^91^92(1 + ^1^2) 

where Vg^g^ is the combined dilution-weighted tag, and Cg-^g^ is the combined efficiency. In 
this manner, tags in agreement as well as tags in conflict are handled properly: in the cases 
where the charge of the two tags agree, the effective dilution is increased; in the cases where 
the two tags disagree, the effective dilution is decreased. 

To help understand the expression for combined dilution V, we examine several limiting 
cases. In the case of a perfect first tagging algorithm, = 1, the combined tag always 
equals the value of the perfect algorithm {Vg^g^ = Vi), independently of the second tagging 
algorithm. For the case where the first tagging algorithm is random, \T>i\ — 0, the combined 
tag always equals the value of second algorithm {Vg^g,^ = 1)2). In the case where the result 
of first tagging algorithm is equal and opposite to the result of the second tagging algorithm 
{Vi = -V2), the V 

qiq2 — 0- This is expected when the two tagging algorithms have equal 
power but give the opposite answer. 

To understand the combined efficiency eg^g^, we consider an example. There are nine 
possible efficiencies for the combined tagging algorithms, Sg^g^. The individual efficiencies 
for perfectly efficient symmetric tagging algorithms have the values e+ = e_ = 0.5 and cq = 
(e_|_-|-e_-|-eo = !)• In this case, five of the nine combined efficiencies are trivially zero. For the 
case of two perfect tagging algorithms giving the opposite result {Vi — —T>2 and [Dil = 1) 
then the combined efficiency must be Cg^g^ — 0, independent of the magnitude of e^^ and 
. This is expected because, by definition, perfect tagging algorithms can not disagree. 
There are only two remaining nonzero cases to examine for the perfectly efficient tagging 
algorithm. For the case in which they agree, the combined efficiencies are e+i^+i = 0.5 and 
e_i,_i = 0.5. 



V. THE LIKELIHOOD FUNCTION 

An extended log-likelihood method is used to determine the best value for sin 2/3, a 
free parameter in the fit. It is helpful to refer to the parameters collectively as a vector 
p with 65 components. The remaining 64 parameters describe other features of the data 
(signal and background) which need to be determined simultaneously, but have only technical 
importance. 

The main ingredient of the likelihood function is the product YliVi where i runs over 
all the selected events and Vi is the probability distribution in the measured quantities: 
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the normalized mass, the flight-time, and the tags (?2, ^s)- The tags, although discrete 
variables, are conceptually thought of as analogous to continuous variables, such as the 
measured mass. The parameters p control the shape of the Pj. There is a separate set 
of parameters for the SVX sample and the non-SVX sample to control the shape of the 
components of Vi. This is especially important for the parts of the function that specify the 
distribution of the measured flight-time and mass, but also the distribution of SST tags. 

The form for Vi assumes that all events are of three types: signal, prompt background, 
and long-lived background. Each possibility is included in Vi. Because the distributions in 
mass, flight-time, and tag are different for the three types, Vi contains separate components 
Vs, Vp, and Vl, which are the overall distributions for signal, prompt background, and long- 
lived background respectively. Additional parameters — a separate set of parameters for SVX 
and non-SVX — specify the relative quantities of each event-type. Each of the components 
Vs, Vp, and Vl is expressed as the product of a time-function (Tg, Tp, Tl), a mass-function 
(Mg, Mp, Ml), and a tagging-efficiency-f unction (£^g, Sp, Si^). 

The time-function Tg is the probability distribution for the observed-time given the 
observed tags, and therefore has a dependence on the measured time and its uncertainty, the 
measured tags and dilutions, and sin 2(3. The lifetime r and mixing parameter Am^ are 
constrained at the world averages: r = (1.54 ±0.04) ps and Am^ = (0.464 ±0.018) hps~^ [Q. 
The Tp function is a simple Gaussian representing the prompt J/i(j background, and depends 
on the measured time and uncertainty. There are two time- uncertainty scale factors in p, 
one for SVX events and one for the non-SVX events, to allow for the possibility that the 
measured time-uncertainties are different from the true uncertainties by a constant factor. 
The Tl function has positive and negative exponentials in time to represent positive and 
negative long-lived background. The positive long-lived background arises primarily from 
real B decays, while the negative long-lived background is used to describe non-Gaussian 
tails in the lifetime resolution. 

The mass-function Mg is a Gaussian representing the normalized mass, and also includes 
a mass-uncertainty scale parameter. The mass-functions Mp and Ml are linear in mass and 
normalized over the ±20cr mass window. 

The tagging-efficiency-function Ss gives the probability of obtaining the observed com- 
bination of tags for a signal event. In addition to the observed tags for the event, it also 
depends on the individual tagging efficiencies and dilutions. The prompt and long-lived 
background tagging-efficiency- functions, Sp and Si,, give the probability of obtaining the 
observed combination of tags for prompt and long-lived background events; they depend 
on individual background tagging efficiencies, but no dilutions are involved because there is 
no right or wrong sign in tagging background. For each individual tagging algorithm, the 
efficiencies and the dilutions (each a component of p) float and are allowed to be different for 
+ and — tags and the corresponding efficiencies and the dilutions for the tag-0 cases follow 
by normalization. However, for the signal, there are constraints on the individual tagging 
efficiencies and dilutions based on the available measurements and their uncertainties. 

A. The likelihood function definition 

The negative log-likelihood i{p) is given by 



15 



B 

1 fm-ifjY 



2 V (Ti 

The 4 free parameters iV|^^, A'^l^'^, jv^on-svx^ ^^^^ jynon-svx j-q^qj. ^j^g number of signal and 
background events in the SVX and non-SVX respectively. The summation over j represents 
a summation over all of the constraints we place on the parameters. The constraints in 
general connect some function fjlp) of the parameters with the corresponding value (fj) 
and uncertainty aj determined by other measurements. 

The summation over i above runs over all data events that satisfy our selection criteria; 
Vi is the probability for the ith event, and implicitly depends on p. The function Vi is given 
by 

= NsVs + Nb [(1 - Fl)Pp + FlPl] 

All events are classified as either type SVX or type non-SVX: the Ns, Nb, and Fl in the 
expression above are actually parameters N^^"^, N^^-^, and F^vx ^^^^q long-lived fraction 
of SVX background) for SVX-type events and AT^on-svx^ jynon-svx^ ^non-svx 
SVX-type events. Although the lifetime resolution for non-SVX events is poor relative to 
the SVX events, the information is used in the likelihood function. 

The functions Vs, Vp, and Vl are the probabilities for the signal, prompt background, 
and long-lived backgrounds. They are given by the products of time, mass, and tagging- 
efficiency functions: 

Vs = TsMsSs Vp = TpMpSp Vl = T^M^^Si, 
The signal time function is specified by 



h{t) = £9192(1 + AcpVg^q^ sin(Amt)) 

r 

where g*h{t) represents the convolution of h{t) with a Gaussian of width a and depends 
implicitly on the values of the flight-time-uncertainty a and sin 2/5. The St above is Sf^-^ 
(the SVX lifetime error scale) for SVX events and 5*^°^^"^^-^ for non-SVX events. The at is 
the uncertainty on the flight-time t of the 5-candidate, determined independently for each 
event. The prompt background allows the determination of Sf^^ and 5*°°°"^^^ using the 
global fit. Knowledge of the individual tag dilutions is incorporated through the constraints. 
The signal mass function is 

where Mb is the normalized mass of the -B-candidate and Sm is the S-mass error scale. 

In an analogous fashion to the combined signal tagging-efficiency function £^S) calcu- 
lated by combining 3 tags as in section |V^, depends on the 8 tagging dilution components 
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(as in Table p|) of p and the 8 individual + and — tagging-efficiency components. The com- 
bined efficiency Ss is the efficiency for obtaining the particular combination of tags observed 
in the event. 

The prompt background time and mass functions are 



a = S,a, St = Sr^" or ^r""'^'' 

Mp = (1 + CMB)/i2W) W = 20 

where W represents the normalized-mass window-size (±20(t), and is the mass-slope of 
the prompt background. 

The combined prompt-background tagging-efficiency function Sp is given by the product 
of the individual prompt background tagging-efficiencies: £p = life £p where k runs over the 
tags. The individual prompt background tagging-efficiencies are parameterized as 

r e^(l- A^)/2 q'' = -1 
£^= < 1-el. q'' = 

[ e|(l + A|)/2 g'^ = 1 

where g'^ is the tag-result of the kth tagging algorithm, and ep and Ap are components of p 
(specifically A^^^^^^, e^^^--^^^, A^^^"""-^^^, e^^^, Af^, ef^^, and Af^^). The 

parameters are the asymmetries of the kth algorithm in tagging the prompt background. 
The SSTsvx and SSTnon-svx are mutually exclusive — k always runs over 3 tags. 
The long-lived time function Tl is given by 



F_^e^l^' t < 

ZT— 

(1 -F_)^e-*/^+ t > 



where F_ is one of F^vx ^^^^^ ^non-svx^ -g ^svx ^^^^^ ^non-svx^ ^^^^ ^_ -g ^svx 

and r°°"-sv^. 

The long-lived mass and tagging-efficiency functions are 

Ml = (1 + ^^MB)/(2iy) £l = 

k 

( eiil - Ai)/2 q^ = -1 
£l= } 1-el g'^ = 

[etil + Ai)/2 g'= = l 

where the notation is exactly analogous to the Mp and Sp defined above. 

To further illustrate the role of constraint terms in the negative log-likelihood function 
we highlight the dilution constraints. There are two dilution parameters, Z)+ and D_, per 
tagging method, the 8 parameters in i{p) representing the tagging dilutions that float in the 
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fit that locates the minimum of (-{p). The probability Vi of the ith J/ipK^ candidate depends 
on these parameters through Ts and Ss- Each tagging method also has its own calibration 
information derived from other decay modes. For example, the dilutions are constrained 
using results from the J/ipK^ calibration sample. In addition, the and D_ dilutions 
for the SST SVX sample are constrained to the average dilution (Z^avc = 16.6 ± 2.2%) 
obtained after extrapolating the mixing analysis dilution to lower Pt [[TOj^ . The available 
calibration information for each tagging method is represented in by constraint terms. 
These terms cause the function to increase as the dilution parameters wander from 
the values preferred by the calibration. When locating the minimum of i{p) we are then 
simultaneously determining sin 2/5 and the 8 dilution parameters, so that the uncertainty on 
sin 2/3 from the fit includes contributions from all of the calibration uncertainties. 

There are similar constraint terms for the efficiency ratios for each tagging method 
(e-|_/e_). The efficiency ratios e+/e_ for each tag algorithm are constrained using the in- 
clusive B J/ipX sample. We fit the J/ip mass distributions for the number of + and — 
tags. The ratio of the number of + tags to the number of — tags constrains e+/e_. The 
B J/ipX sample is assumed to have negligible intrinsic CP asymmetry. In addition, the 
B^ lifetime tbo and mixing parameter Am^ are free parameters in the fit, and there are terms 
to constrain each to its world average 0. The parameter tbo is constrained to 1.56 ±0.04 ps 
and the parameter is constrained to 0.464 ± 0.018 ^ps~^. Although constraining Am^ 
to the world average is the most natural procedure, we also have the option of determining 
Amrf and sin 2/5 simultaneously from the J/ipKg data by removing the constraint on Am^. 

The calibration measurements are summarized in Table |T|. The efficiency ratios are 
consistent with expectations. For SST, the ratios are greater than unity due to a higher 
efficiency for reconstructing tracks with positive charge in the CTC. 



B. Fits to toy Monte Carlo data 

As a check of the fitting procedure several sets of ~ 1000 toy Monte Carlo data samples 
were generated, each set generated with a different value of sin 2(3. The number of events, 
SVX/non-SVX ratio, signal-to-background ratios, tagging efficiencies and dilutions, mass 
uncertainty and its scale factor, background lifetimes, time uncertainties and scale factors, 
and other kinematic features of the generation procedure were all tuned to be similar to the 
composition of the data sample. 

The left plot in Fig. ^ shows the distribution of the appropriate uncertainty (allowing 
for asymmetric errors p9[) on sin 2/? returned from the Monte Carlo fits with generated 
sin 2/? = 0.5. The typical value of the uncertainty on sin 2/5 returned from these fits is 
~0.44, though there is a long tail extending out to ~0.7. The width of the distribution is 
determined by Poisson fluctuations in the number of Monte Carlo events that are tagged. 
The right plot in Fig. ^ shows [sin2/?(flt) — 0.5]/a, where a is the appropriate + or — 
uncertainty on sin 2/5. 

The results from this and other samples generated at different values of sin 2(3 support 
that the fltting procedure provides an unbiased estimate of the value of sin 2(3 of the parent 
distribution. The distribution of the difference between the flt-sin2/5 and the true sin 2/5 of 
the parent distribution is well approximated by a Gaussian and the flt-uncertainty on sin 2(3 
provides a good estimate of the a of that Gaussian. 
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C. Systematic uncertainties 



Systematic uncertainties on the measurement of sin 2/3 due to flavor tagging, the B 
hfetime and Am^ are included as constraints in the fit. We evaluated the systematic uncer- 
tainties due to the uncertainty in the mass, trigger bias and regeneration. 

The systematic uncertainty arising from the B mass is studied using 1000 simulated 
experiments. The data were generated at the nominal B mass and three full likelihood 
fits were performed on each experiment. One fit was performed using the normalized mass 
calculated with the nominal B mass and two additional fits were performed using B masses 
shifted by ±1 MeV/c^. The shifts observed in sin 2/3 from fits to the simulated experiments 
are consistent with a random distribution centered on zero with an RMS of 0.019. The 
change in the observed RMS spread of sin 2/? is < 0.019 when combined in quadrature. 
We also fit the data with the B mass shifted by 1 MeV/c^ and found the value of sin 2/? 
changed by 0.013, which consistent with the simulation results. We conclude the additional 
uncertainty on sin 2/? due to the uncertainty on the B mass is < 0.019 and is negligible. 

The data are assumed to be a 50:50 mix of B^ / B^ . A possible charge bias arising from 
the trigger is considered. Events that are triggered on the two muons from the J/ip decay 
do not contribute to the charge bias. The remaining 30% contain some events in which the 
trigger was from one of the J / ip muons and the other lepton candidate was from the opposite 
side B. The magnitude of the charge bias in the trigger has been measured to be < 1% at 
a threshold of Pt = 2 GeV/c and is consistent with zero for Pt > 3 GeV/c, rendering this 
uncertainty negligible. 

Possible contamination of our data from regeneration from the material in the inner 
detector has been considered. Reconstruction of the as a Kg causes the event to be 
entered with the incorrect sign in the asymmetry. This effect shifts sin 2/? by less than 
0.003, which is neglected. The results of the systematic studies are shown in Table |IV|. 

We have evaluated the contribution to the sample from B^ — ^ J/ipK*, with K* K^tt^ 
and the 7r° not reconstructed and find it to be a negligible contribution. The same is true 
with — >■ J/ipA and A — pvr" and the A reconstructed as Kg — > tt+tt"; Bs — >■ J/ip(f), 
^ K'^sK'^ and Bs J/^K% 

Many checks of the data and analysis have been performed to increase our confidence 
in the result. In order to check the sensitivity of the result to the dilutions, we imposed 
alternative JETQ and SLT dilution parameters taken from our various mixing analyses that 
use the inclusive lepton sample We observe the expected shift in the value of sin 2/? 
and small changes in the uncertainty. The signal sample selection criteria have been varied, 
and other than a sensitivity to the SST tag track Pt threshold, as discussed in Ref. |]10|, we 
find no unexpected sensitivity in the result. 



D. Final result 

The maximum likelihood function fitting procedure returns a stable value for sin 2/? 
and the uncertainties are approximately Gaussian. Even though asymmetric dilutions are 
permitted in the fit, no significant asymmetry is observed. Furthermore, the background 
asymmetries are consistent with zero. 

Using the entire data set and three tagging algorithms, we find 
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sin 2/5 = 0.79 

^ -0.44 

The asymmetry is shown in Fig. ^ for the SVX and non-SVX events separately. The 
asymmetry for the SVX events is displayed as a function of lifetime, while the asymmetry 
for the non-SVX events is shown in a single, time-integrated bin since the decay length 
information is of low resolution. Although plotted as a time-integrated point, lifetime infor- 
mation for the non-SVX events is utilized in the maximum likelihood function. The positive 
asymmetry preferred by the fit can be seen. The curves displayed in the plot are the results 
from the full maximum likelihood fit using all data. In order to display the data, we have 
combined the effective dilution for single and double-tag events after having subtracted the 
background. The full maximum likelihood fit uses the SVX and non-SVX samples and treats 
properly the decay length, dilution and uncertainty for each event. 

The uncertainty can be divided into statistical and systematic terms: 

sin 2/? = 0.79 ± 0.39(stat) ± 0.16(syst). 

The systematic term predominantly reflects the uncertainty in the result due to the un- 
certainty in the dilution parameters. Although the dilution parameters are not precisely 
determined, due to the limited statistics of the J/ipK^ calibration sample, this 

uncertainty term does not dominate the overall uncertainty on sin 2/5. Furthermore, the 
uncertainty on sin 2/? will not be dominated by the uncertainty on the dilution parame- 
ters in future runs because the uncertainty scales inversely with increasing statistics of the 
calibration samples. 

It is of interest to determine the quantitative statistical significance of whether this result 
supports sin 2/5 > 0.0 and hence provides an indication of CP symmetry violation in the b 
quark system. A scan through the likelihood function as sin 2f3 is varied is shown in Fig. |^ 
and demonstrates that the uncertainties follow Gaussian statistics. Using the Feldman- 
Cousins frequent ist approach |3^, we calculate a confidence interval of 0.0 < sin2/? < 1 at 



93%. An alternative approach is the Bayesian method, where a flat prior distribution in 
sin 2/5 is assumed and a probability that sin 2(3 > 0.0 of 95% is calculated. Finally, if the true 
value of sin 2/5 is zero, and the measurement uncertainty is 0.44 (Gaussian uncertainty), the 
probability of obtaining sin 2/5 > 0.79 is 3.6%. This value is obtained by simply integrating 
the Gaussian distribution from 0.79 to oo. The toy Monte Carlo is in good agreement with 
the calculated probability. 

It is possible to remove the constraint that ties Am^ to the world average value and 
to fit for sin 2/5 and Am^ simultaneously. In this case the result is sin 2/5 = 0.88+J;JJ and 
Arud = 0.68 ± 0.17 ^ps~^. The value of Am^ from the fit agrees with the world value at 
the level of ~ 1.2a. This agreement increases our confidence in the main result. Figure |^ 
shows the la "error ellipse" contour in sin 2/5- Am^^ parameter space for the fit when both 
parameters float freely, and for comparison the nominal sin 2/5 result with the world average 
Arrid and uncertainty. From the roughly circular shape of the contour, the Am^ and sin 2/5 
parameters are largely uncorrelated in the fit. 

A time-integrated measurement to check the final result was performed. This simplified 
analysis does not use the time dependence of the asymmetry and ignores the small tagging 
asymmetry corrections applied in the full maximum likelihood fit. Each event falls into 
one of 12 classifications depending upon the type of flavor tags available for that event. 
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Each event can be associated with only one class of tag combination. The effective tagging 
efficiency for the entire sample, eD^, is (6.3 ± 1.7)%. A value of sin 2/3 for each class is 
calculated and a weighted average from the 12 classes is determined. Ignoring correlations 
in the dilution, sin 2/5 = 0.71 ± 0.63. This value is consistent with the final result and 
demonstrates the improvement in the uncertainty of sin 2/3 provided by the full maximum 
likelihood procedure. This improvement agrees well with improvements observed using the 
toy Monte Carlo. 

Table summarizes fit results for various tag-dataset combinations. The three tagging 
algorithms contribute roughly equally to the precision of the sin 2/3 measurement. Although 
the SVX and non-SVX sample sizes are approximately equal, the SVX events contribute 
more significantly to the final result. The main reasons for this are that the precision 
lifetime information from the SVX allows a better determination of where the decay takes 
place along the oscillation curve and the better signal-to-background level from eliminating 
the prompt background. 

The row in Table labelled SVX SST is the result obtained when this analysis restricts 
the data set to the SVX sample and uses only the SST algorithm. This procedure essentially 
repeats the published CDF sin 2/3 analysis that obtained sin 2/3 = 1.8 ± l.l(stat) ±0.3(syst). 
The small difference is due to sample selection. 

VI. MIXING IN THE B J/ipK* SAMPLE AS A CHECK 

A control sample of 5° J/tljK*{892f decays, where K*{892)^ K^n^, can be 
analyzed for the presence of an oscillation due to mixing (Am^ is well measured) in order 
to verify the tag algorithms and likelihood fitting procedure. The three flavor tagging 
algorithms are used to determine the neutral B flavor at the time of production and the 
dilution parameters are constrained using the same values as in the B J/ipKg analysis. 
The charge of the kaon is used to differentiate the 5" from b'^ at the time of decay. After 
correcting for tagging dilutions, the amplitude of the oscillation still differs from unity due 
to the probability that the K^tc' is reconstructed as K~tt'^, which occurs about = 5% 
of the time due to the wide K* resonance. 

The J / ip K g- J /ipK* (892)^ analogy is however not perfect. In order to achieve similar 
signal-to-background ratios, the selection criteria for the B J /ipK* (892)^ are more severe, 
which changes the kinematic properties of one sample with respect to the other. The largest 
backgrounds for both decay modes are at short decay distances and they decrease as the flight 
path increases. This works to our advantage in the CP analysis but reduces the sensitivity 
of the mixing analysis. In particular, due to the different oscillation phase in the CP analysis 
versus this mixing analysis, (sin(Amdt) —>■ coslArridt)) , the smallest signal-to-background 
ratio occurs at the peak of the mixing amplitude for B J /ipK* {892Y data set, where as a 
very favorable signal-to-background ratio occurs at the peak of the B J / ipK^ oscillation. 
In both the J/ipKg and J/ipK* (892)^ modes, 75-80% of the background is prompt, i.e. 
consistent with having zero lifetime. 

The sample is constructed using similar criteria to that used to reconstruct the B — >■ 
J/ipK decay modes in this paper. The J/ip selection for this decay mode is the same 
as the J/ipKg analysis. Pion and kaon tracks are required to have Pt > 500 MeV/c. The 
reconstructed i^*(892)° candidates are required to have an invariant mass within 80 MeV/c^ 
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of the world average of 896.10 ± 0.28MeV/c2 H i^*(892)° mass. The K* candidate must 
have Pt > 3 GeV/c. The four-track fit for J /ifjK* is the same as the fit for J/ipK^, except 
the four tracks are required to meet at a common point and the K* mass is not constrained. 
If a candidate event has two tracks that satisfy two K*{892)^ combinations {K^tt~ / K~n^) 
then the combination with a Kit mass closest to the mean i^'*(892)° mass is chosen. Finally, 
if multiple K* candidates are found in an event, the /('*(892)° candidate chosen is the one 
that gives the best four-track fit. All four charged tracks {fi,fi,K,n) must originate from a 
common vertex and a Pt{B) > 4.5 GeV/c is required. A total signal sample of 226 ± 24 
events where both muon candidates have precision lifetime information and 231 ± 28 events 
where < 1 muon candidate has precision lifetime information are found. 

The maximum likelihood fit to the J/iIjK*{892)^ data is implemented in the same way 
as previously described for J/ipKg except for the time- function Ts in which h{t) is replaced 
by: 

Q-t/r 

h{t) = 69192(1 + T^kT^qiqi cos(Amdt)). 

r 

Here Vk = QkOki where qk is the charge of the K"^ from the decay of the i^*(892)°, 
and Dk is the dilution arising from the inability to correctly distinguish the charged kaon 
from the charged pion in the ii'*(892)° decay. The dilution is the free parameter in this 
fit and is analogous to sin 2/3 in the J/ipK^ fit, the parameters in each case representing the 
amplitude of an oscillation. The amplitude is expected to be Dk = 1 — 2Pk = 0.9^^^ where 
Dk is the dilution factor coming from incorrect K-ir assignment 



When Arud is fixed to the world average, we measure Dk = 1.00 ± 0.37, which is 
consistent with expectation. When Ama is allowed to float, we measure: Dk = 0.96 ± 0.38 
and Arrid = 0.40 ± 0.18h ps~^, which is consistent with the world average Arrid = (0.464 ± 
0.018) hps-^ II . The results of the fits are shown in Fig ^ 

Although the statistics are not sufficient for a precise measurement of Am^, this check 
on an independent sample of events is entirely consistent with our expectation. 



VII. CONCLUSION 

We have presented a measurement of sin 2/3 using ~ 400 B J/tpKg events recon- 
structed with the CDF detector. We find: 

sin 2/5 = 0.79'^U'^^(stat + syst) 

with the uncertainty dominated by the statistical contribution. 

We have calculated the statistical significance of whether this result supports sin 2(3 > 0.0 
and hence provides indication for CP symmetry violation in the b quark system. Using the 
Feldman- Cousins ||3^ method, a 93% confidence interval of 0.0 < sin 2/3 < 1.00 is found. 
Alternative methods yield similar limits. This measurement is the best direct indication that 
CP invariance is violated in the b quark system and is consistent with the Standard Model 
expectation of a large positive value of sin 2/3 IH^-Q. The sign of our result supports the 
favored positive signs for Bb and Bk- With an anticipated luminosity of 2fb~^ in Run II, 
we expect, based on a simple extrapolation of this measurement, an uncertainty on sin 2/? 
of ~0.08. Detector upgrades in progress should further reduce this uncertainty. 
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FIG. 3. Left: Normalized mass distribution of the J/ipKg candidates where both muons have 
good SVX information providing a high precision decay length measurement. Right: Normalized 
mass distribution of the J/tpK^ candidates in the non-SVX sample. Either one or both muons 
are missing good SVX information, leading to a low resolution decay length. For both plots, the 
curves are Gaussian signals plus linear background. 
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FIG. 4. The mass distribution of the J/ijjK"^ candidates both with and without SVX informa- 
tion. The curve is a Gaussian signal plus linear background from the likelihood fit. 
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FIG. 5. Left: Distribution of C7sin2^ from fits to multiple Monte Carlo datasets generated with 
sin 2/3 = 0.5. Right: Distribution of normalized sin 2/3 deviations, i.e. (fit-sin 2/3 — 0.5)/(7sin2^, and 
a Gaussian fit to that distribution. The mean of the Gaussian fit is 0.038 it 0.033 and the width is 
1.01 ± 0.03, consistent with expectation. 
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FIG. 6. The true asymmetry (sin 2/3 sin Am^t) as a function of lifetime for B — > J/tjjKg events. 
The data points are sideband-subtracted and have been combined according to the effective dilution 
for single and double-tags. The non-SVX events are shown on the right. 
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FIG. 7. A scan of the log-likelihood function. The value of sin 2/3 is scanned, and at each step, 
the function is minimized. 
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FIG. 8. The la (39%) sin2/3-Amrf contour from a fit with ^rrid constrained only by the 
J/'4)K'^s data. Also shown is the nominal fit with Am^ = (0.464 ± 0.018) /ips~^ g. 



33 




0.05 0.1 0.15 0.2 0.25 

ct (cm) 

FIG. 9. The true asymmetry {Dk cos Amat) as a function of lifetime for 5° J/il;K*{892)^ 
events. The data points are sideband-subtracted and have been combined according to the effective 
dilution for single and double-tags. The time-integrated asymmetry for non-SVX events is shown 
on the right. The solid curve represents the maximum likelihood fit in which Arrid is fixed and the 
dashed curve is the expectation when we also fix Dk- 
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TABLES 



TABLE I. Summary of tagging algorithms performance. All numbers listed are in percent. 
The efficiencies are obtained from the B — ^ J/tpKg sample. The dilution information is derived 
from the — > J/ijjK^ sample. 



tag side 


tag type 


class 


efficiency 


dilution 


same-side 


SST 


in SVX 


35.5 ± 3.7 


16.6 ±2.2 




SST 


/il or /i2 non-SVX 


38.1 ±3.9 


17.4 ±3.6 


opposite side 


SLT 


all events 


5.6 ± 1.8 


62.5 ± 14.6 




JETQ 


all events 


40.2 ± 3.9 


23.5 ± 6.9 



TABLE IL Definition of tags. For the case of the SST algorithm, the tag depends upon the 

charge of a track (t^ ,t~) near the B; for the SLT algorithm, the tag depends upon the charge of a 
lepton in the event for the JETQ algorithm, the tag depends upon the average weighted 

charge of tracks in a jet (Qjct)- 



tag positive (+) tag 


negative (— ) tag 


no tag 


SST single track t+ 
SLT single lepton £~ 
JETQ Qjet < -0.20 


single track t~ 
single lepton .£+ 
Qjet > 0.20 


no track 
no lepton 

IQjetl < 0.20 


TABLE in. The dilutions determined from the B^ J/ipK^ sample and 
determined from the inclusive J/'0 sample are shown. Dave is the average 
dilutions utilize additional information as described in the text. 


the efficiency ratios 
dilution. The SST 


tag 


D+(%) D-{%) 


i?ave(%) 


SSTsvx 1.031 ± 0.011 
SSTnon-SVX 1.037 ±0.010 
SLT 0.978 ± 0.047 
JETC) 0.977 ± O.OL") 


16.1 ±5.1 17.1 ±5.2 
17.0 ±5.7 17.8 ±5.8 
76.9 ± 19.6 46.4 ±21.8 
20.7 = 9,:-! 20.5 = S.:-! 


16.6 ±2.2 

17.4 ± 3.6 
62.5 ± 14.6 
2:i.5 = ().9 


TABLE IV. Systematic uncertainties in the measurement of sin 2/3. The items labelled "in 
fit" are parameters that are allowed to float in the fit but are constrained by their measured 
uncertainties. The uncertainty returned from the likelihood fit includes the contributions from 
these sources. 


parameter 


(5 sin 2(5 


in fit 


dilution and efficiency 
Anid 

niB 
trigger bias 
regeneration 


0.16 
negligible 
negligible 

negligible 
negligible 
negligible 


yes 
yes 
yes 
yes 
no 
no 



35 



TABLE V. Fit sin 2/3 results for the three tagging algorithms. The combined for the SST, 
JETQ, and SLT tagging algorithms is 4.63 for 2 degrees of freedom, giving a probability of ~ 10%. 



data tag(s) sin 2/3 + error — error 

ail all 079 041 044 

SST 2.03 0.84 0.77 

JETQ -0.31 0.81 0.85 

SLT 0.52 0.61 0.75 

SVX all 0.54 0.52 0.57 

SST 1.77 1.04 1.01 

non-SVX all 1.24 0.75 0.70 
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